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Abstract 

We propose cosmological redshift-space distortion of correlation functions of galaxies and 
quasars as a probe of both the density parameter VLq and the cosmological constant Aq. 
In particular, we show that redshift-space distortion of quasar correlation functions at 
z ^ 2 can in principle set a constraint on the value of Aq. This is in contrast This is in 
contrast to the popular analysis of galaxy correlation functions in redshift space which 
basically determines fig^/fe, where h is the bias parameter, but is insensitive to Aq. For 
specific applications, we present redshift-space distortion of correlation functions both in 
cold dark matter models and in power-law correlation function models, and discuss the 
extent to which one can discriminate between the different Aq models. 



Subject headings: cosmology: theory — large-scale structure of the universe — methods: 
statistical 



1 Introduction 



Redshift-space distortion of galaxy two-point correlation functions is known as a power- 
ful tool in estimating the cosmological density parameter Qq; on nonlinear scales Davis 
& Peebles (1983) computed the relative peculiar velocity dispersions of pair of galaxies 
around lh~^Mpc from the anisotropies in the correlation functions of the CfAl galaxy 
redshift survey, and then concluded that Qq = 0.20e^°'^ (also see Mo, Jing & Borner 
1993; Suto 1993; Ratcliffe et al. 1996). In linear theory Kaiser (1987) showed that the 
peculiar velocity field systematically distorts the correlation function observed in redshift 
space; the averaged redshift-space correlation function ^^^\x) of galaxies is related to its 
real-space counter part ^^^\x) as 

e(^)(x) = fl + |/5 + Me('^)(a;), (1) 
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where b is the bias parameter, f{z = 0) is the logarithmic derivative of the linear growth 
rate with respect to the scale factor a (eq. |TB[ below) at z = 0. As is clear from the 
empirical fitting formula in the second line by Lahav et al. (1991), however, this formula 
depends on but is practically insensitive to the cosmological constant Aq. 

At higher redshift z ^ 1, another anisotropy due to the geometrical effect of the spa- 
tial curvature becomes important. We derive a formula for the cosmological redshift-space 
distortion in linear theory. This formula turns out to be a straightforward generalization 
of that derived by Hamilton (1992) for z = 0. It provides a promising method to estimate 
both Qq and Aq from future redshift surveys of galaxies and quasars including the Sloan 
Digital Sky Survey (SDSS) and the 2dF. A very similar idea was put forward indepen- 
dently by Ballinger, Peacock and Heavens (1996) in which they considered the anisotropy 
of the power spectrum while we developed a formulation in terms of two-point correlation 
functions. We outline the derivation in §2, and present some examples of results in cold 
dark matter (CDM) models and in power-law correlation function models. 

2 Cosmological redshift distortion 

Throughout the present analysis, we assume a standard Robertson - Walker metric of the 
form: 

ds^ = -dt^ + a{ty{dx^ + Six^ide^ + sin' Od(j)^]}, (3) 
where S{x) is determined by the sign of the spatial curvature K as 

sm{VKx)/VK {K>0) 
Six) = { X iK = 0) (4) 

smh{y/^x)/V^ {K<0) 
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In our notation, K is written in terms of the scale factor at present oq and the Hubble 
constant Hq as 

K^alHl{^l^ + Xo-l). (5) 
The radial distance x{^) is computed by 

where we introduce the Hubble parameter at redshift z: 



H{z) = H^^Jn^il + zf + (1 - Qo - Ao)(l + zf + Ao. (7) 

Consider a pair of galaxies or quasars located at redshifts zi and Z2. If both the 
redshift difference 8z = zi — Z2 and the angular separation of the pair 59 are much less 
than unity, the comoving separations of the pair parallel and perpendicular to the line- 
of-sight direction, and x^, are given by 

where c^\{z) = Hq/H{z), c±{z) = aoHQS{x{z))/z and z = {zi + Z2)/2. Thus their ratio 
becomes 

x^{z) c^{z) z5e ~ ' z5e' ^ ' 

Since X\Jxx_ should approach 5z/{z59) for 2; <^ 1, ri{z) can be regarded to represent the 
correction factor for the cosmological redshift distortion (upper panel of Fig.l). This was 
first pointed out in an exphcit form by Alcock & Paczyhski (1979). 

Although r){z) is a potentially sensitive probe of Jig and especially Aq, this is not 
directly observable unless one has an independent estimate of the ratio X\Jx±. Therefore 
Alcock & Paczynski (1979) proposed to apply the result to intrinsically spherical structure 
resulting from gravitational clustering. More recently Ryden (1995) proposed to apply 
the test to statistically spherical voids around z ^ 0.5. It is not likely, however, that the 
intrinsic shape of the objects formed via gravitational clustering is spherical even in a 
statistical sense. Thus their proposal would be seriously contaminated by the aspherical 
shape distribution of the objects considered. In what follows, we propose to use the 
clustering of quasars and galaxies in z, which should be safely assumed to be statistically 
isotropic, in estimating r]{z). Also we take account of the distortion due to the peculiar 
velocity field using linear theory. 

If an object located at redshift Zi has a peculiar recession velocity Vy^i with respect to 
us, we actually observe its redshift as 

1 + ^obs,l = (1 + ^1) (1 + t^ii,! - 1^11,0) (10) 
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where f y^o is our (observer's) peculiar velocity, and we assume that all peculiar velocities 
are non-relativistic. Let us consider a sample of galaxies or quasars located in the range 
z — dz ^ z + dz with dz z distribute over an angular radius comparable to the proper 
distance corresponding to dz. Then we set the locally Euclidean coordinates in real 
space (comoving), x, and in redshift space, s, which share the same origin at redshift z. 
Since we assume that these coordinates are far from us, we can use the distant-observer 
approximation. If we take the third axis along the line of sight, the above two coordinates 
are related to each other as 



S3 



Xi 



c±[z) 

ZohsA — Z 



S2 



X2 

c±{z 

1 







Cii U 



l + z. 



^011 



(11) 

(12) 



Note that in the last expression, we assume that X3 is sufficiently small and replace 
1 + z + Hxs by 1 + z. The number densities of objects in these coordinates are related by 
the Jacobian of the above transformation. In linear order of density contrast and peculiar 
velocity, we obtain 

5^^\six)) = 5^^\x)-^dsv,. (13) 

According to linear theory, the peculiar velocity is related to the mass density fluctuation 
6m (e.g., Peebles 1993) as 



VAX 



H{z) 
l + z 



f{z)d,A-X,{x) 



where A is the inverse Laplacian, 



d\nD{z) 
din a 



Qiz\ 



,0.6 



+ 
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Q{z) 



and the linear growth rate D{z) (normalized as D{z) = 1/(1 

l + z' 



z) for z — > 00) is 



D{z) = ^^H{z) 



H(z'Y 



dz'. 



(14) 



(15) 



(16) 



In the above, Q{z) and X{z) are the density parameter and the dimensionless cosmological 
constant at z: 



n{z) 



Ho 
H(z) 



[l + zfQo, A(z) 



Hp 
Hiz) 



Ao. 



(17) 



Allowing for the possibility that the density contrast of mass, 6^, differs from that of 
objects 5^*"^ we assume linear biasing S^'''^ = b Sm, where b{z) is a bias factor which now 
depends on z. Then equation (|13[) simply generalizes Kaiser (1987) 's result to z 7^ 0: 



6^'\s{x)) 



d^k 
(2^ 



1 + /^W^ 



ikx 
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where P{z) = f{z)/b{z), and (5^^-' is the Fourier transform of the density contrast. From 
equation (|18D, we derive the following relation for the two-point correlation function at 
2; 7^ which also generalizes the formula of Hamilton (1992) at z = 0: 

^^'\s±, s,) = (1 + + lim?) ux)Po{f^) - [Im + ^im?) 6(x)P2(/i) 

+ hm]%ix)P,i^^), (19) 

60 

where x = \Jc\\^S\^'^ + c^^sj_2, /i = c^\S\Jx (sy = S3 and s]_ = sl + s^), -Pn's are the Legendre 
polynomials, and 

1 f°° 

^2^^) = ^ dkk'^j2i{kx)P{k;z) 

Again in linear theory and for linear biasing, the power spectrum P{k] z) of objects at z 
is related to that of mass P^^\k) of matter at 2; = as 

^) = (§|y) P^-'Hk), (21) 

3 Simple model predictions 

The lower panel in Figure 1 plots f{z). Here and in what follows, we use the fitting 
formula (|l^) for f{z) whose mean error is 2 percent. As expected from equation (H), 
f{z = 0) is insensitive to the value of Aq and basically determined by Qq only. At higher 
redshifts, however, f{z) becomes sensitive to Aq as pointed out earlier by Lahav et al 
(1991), especially if Qq -C 1. Therefore the behavior of f{z) at low and high z is a 
potentially good discriminator of Qq and Aq, respectively. What we propose here is that 
a careful analysis of the redshift distortion in correlation functions of galaxies at low 
redshifts and quasars at high redshifts can probe Qq and Aq through f{z) as well as ri{z). 

For specific examples, we compute ^^'^\s±, Sy) in linear theory applying equations ([T9|) 
and ( PU[ ) in CDM models with Hq = 70 km-s~^-Mpc~^ (Tanvir et al. 1995). The resulting 
contours are plotted in Figure 2. The four sets of values of Qq and Aq are the same as 
adopted in Figure 1, and we adopt the COBE normalization (Sugiyama 1995; White and 
Scott 1995). 

In order to quantify the cosmological redshift distortion in Figure 2, let us introduce 
the anisotropy parameter ^^f\s)/C,±\s), where ^±\s) = ^''^\s,0) and ^|f''(s) = ^'^''^(0, s). 
The left four panels in Figure 3 show the anisotropy parameter against z in CDM models. 
Since the evolution of bias is largely unknown, we assume 6 = 1 and 2 just for definiteness. 
The effect of the evolution of bias (e.g., Fry 1996) will be considered elsewhere (Suto & 
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Matsubara 1996). This clearly exhibits the extent to which one can discriminate the 
different Aq models on the basis of the anisotropics in ^^^^ at high redshifts. 

For comparison, let us consider a simple power-law model ^'•'"•'(x; z) = A{z)x~'^ (7 > 0). 
Then equation (|19D reduces to 



^^"H\AMTc5i;2) 3-7 (3-7)(5-7) 

In this case the anisotropy parameter is given by 

eiP(^) ^ . , (3 - 7)(5 - 7) + 2(1 - 7)(5 - 7)/3(^) + (3 - 7)(1 - l)^' 

' (3-7)(5-7) + 2(5-7)/3W+3/3(;^)2 ' ^ ' 

and is independent of the scale s. 

The right four panels in Figure 3 show that the behavior of the anisotropy parameter 
is sensitive to the power-law index 7. This is partly because fio = 1 CDM models 
which we consider are already nonlinear on (10 ~ 20)/i~^Mpc (Fig. 2), and then the linear 
theory prediction is not reliable enough. On the other hand, 17 = 0.1 CDM models are 
well described in linear theory on the scales of our interest ^ l/i~^Mpc. Since most 
observations do suggest that Qq in our universe is around (0.1 ~ 0.3) (e.g., Peebles 1993; 
Suto 1993; Ratcliffe et al. 1996), this is encouraging for our present result on the basis of 
linear theory. Incidentally, Figure 3 also implies that = 0.1 CDM models in the linear 
regime are well approximated by the power-law model with 7 ~ 1. This is reasonable 
since the CDM spectral index of P{k) is around —2 which corresponds to 7 ~ 1. 

4 Conclusions 

Redshift-space distortion of galaxy correlation functions has attracted much attention as 
a tool to determine Q^'^/b (Kaiser 1987; Hamilton 1992). We obtained an expression 
to describe the cosmological redshift-space distortion at high z in linear theory taking 
proper account of the spatial curvature K. Then we showed that in principle this can 
discriminate the value of Aq via the 2r-dependence of the l3{z) and r]{z) parameters. 

Recent analysis, for instance of the Durham/UKST galaxy redshift survey on the basis 
of eq.(0), yields n^^Vb = 0.55 ± 0.12 (Ratcliffe et al. 1996) with ~ 2500 galaxies. Their 
Figure 4(a) clearly exhibits that the observational data at z = are already statistically 
reliable for the direct quantitative comparison with the upper panels in our Figure 2, 
although we do not attempt it at this point. SDSS, for instance, is expected to have a 
million of galaxy redshifts up to ^ ^ 0.2 which would be able to determine i^Q^/fo to better 
than 10 percent accuracy. Again upon completion of SDSS, ~ 10^ quasar samples become 
available and the anisotropy in quasar correlation functions at z ~ 2 will put a constraint 
on Aq, given Qq determined from the galaxy correlation functions. The quasar luminosity 
function of Boyle, Shanks & Peterson (1988) predicts that the number of quasars per 
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TT steradian brighter than 19 B magnitude is about 4500 either in 2; = 0.9 ~ 1.1 or in 
^ = 1.9 ~ 2.1 (for the Qq = 1 and Aq = model). Therefore we expect that the resulting 
statistics is even better than that obtained by Ratcliffe et al. (1996) if unknown evolution 
effects and other potential systematics interfere. 

Although the linear theory which we used throughout the paper becomes less restric- 
tive at higher redshifts, it is observationally easier to detect clustering features in the 
nonlinear regime. Thus the analysis of nonlinear redshift-space distortion is another im- 
portant area for research (Suto & Suginohara 1991; Matsubara & Suto 1994; Suto & 
Matsubara 1994; Matsubara 1994). In addition, it is important to examine the possi- 
ble systematic errors due to the finite volume size and the shape of the survey region. 
This can be best investigated by the direct analysis of the numerical simulations. This is 
partly considered by BaUinger, Peacock, & Heavens (1996) although in k-space. We plan 
to return to this in the later paper (Magira, Matsubara, & Suto 1996). 

We thank Naoshi Sugiyama for discussions. This research was supported in part by 
the Grants-in-Aid by the Ministry of Education, Science and Culture of Japan (07740183, 
07CE2002). After we submitted the present paper, we realized the similar independent 
work by Ballinger, Peacock, & Heavens (1996). We are grateful to John Peacock for 
calling our attention to this paper and for useful comments. 
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Figure 1 : r){z) — C\i{z)/c±{z) {upper panel), and f{z) — P{z)b{z) {lower panel) for 
(Qq, Aq) = (1.0,0.0) in dashed line, (1.0,0.9) in dot-dashed line, (0.1,0.9) in dotted 
line, and (1.0,0.9) in thick solid line. 

Figure 2 : The contours of C*'''''("S±, Sy) in CDM models ai z — {upper panels) and 
z — 3 {lower panels) . Hq — TOkm-s^^-Mpc"^ is assumed, and the amplitude of the 
fluctuation spectrum is normalized according to the COBE 2 yr data. Solid and 
dashed hues indicate the positive and negative ^^^\ respectively. Contour spacings 
are Alog^olCl = 0.25. Thick lines in all VLq = 1 models represent ^'^^'^ = 10, 1, 0.1, 
-0.1, and -1. Thick hues in Qq = 0.1 models indicate that = 0.01 and 0.001 
for Ao = at ^ = 0, = 1, 0.1 and 0.01 for Aq = 0.9 at z = 0, = 0.01 for 
Ao = at ^ = 3, and = 0.1 for Aq = 0.9 at ^ = 3. 

Figure 3 : The anisotropy parameter ^\\^\s) / ^^f^ {s) as a function of z. Upper and 
lower panels assume that b — 1 and 2, respectively. Prom left to right, the panel 
corresponds to CDM at s = 10/i~^Mpc, CDM at s = 20/i~^Mpc, a power-law 
model with 7 = 1.8 and a power-law model with 7 = 1.0. 
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fio = 1.0, Xq = 0.0 Qq = 1.0, Xq = 0.9 Qq = 0.1, Xq = 0.0 = 0.1, Xq = 0.9 
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